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ABSTRACT 



This thesis computerizes aircraft Naval Air Training and 
Operating Procedures Standardization Program (NATOPS) flight 
performance charts and shows it to be feasible with a high 
degree of accuracy. The computer programs developed are 
Normal Take-off and Cruise Performance for the A-6E aircraft 
and are adaptable to hand-held calculators, desk calculators, 
or existing aircraft systems computers. The feasibility, 
procedures, and techniques shown in this thesis are applicable 
to any aircraft, both fixed and rotary wing. 

A much improved and more accurate method of mission plan- 
ning is needed to reduce NATOPS chart errors during pre-flight 
planning and to reduce the accident potential of air crews who 
make performance decisions during flight from past experience 
or old information, without accounting for all variables which 
affect aircraft performance. A large percentage of Navy/ 

Marine Corps aircraft accidents has been directly or indirectly 
attributed to misinterpretation of NATOPS charts and/or poor 
performance decisions. Computerized NATOPS performance charts 
will significantly reduce human errors inherent to visual 
chart interpretation and, when coupled with aircrews having 
computer access to make airborne aircraft performance deci- 
sions based on all performance variables, will greatly enhance 
safety of flight in completing any or all mission phases, 
thereby increasing operational readiness of all Navy/Marine 
Corps fleet units. 
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I. INTRODUCTION 



The objective o£ this thesis was to investigate the feas- 
ibility of computerizing aircraft Naval Air Training and 
Operating Procedures Standardization Program (NATOPS) flight 
performance charts with a high degree of accuracy and to 
develop computer programs which are adaptable to hand-held 
calculators, desk calculators, or existing internal aircraft 
systems computers. The NATOPS performance charts of opera- 
tional aircraft used in the Navy and Marine Corps are a 
conglomeration of graphs, charts, curves, and tables which 
are frequently cumbersome and very susceptible to user error 
in their interpretation and manipulation. The charts present 
performance information for all aspects of flight operations, 
from take-off to final landing. Therefore, it is imperative 
that flight crews be able to interpret and utilize these 
charts to obtain optimum and efficient performance from their 
respective aircraft in all phases of the specified mission. 

Mission planning is an essential and key element of flight 
operations and must insure a high degree of safety and maximum 
utilization of the aircraft performance in order to success- 
fully accomplish the mission. This planning requires use of 
the NATOPS charts to determine take-off distances, climb 
profiles, cruise distances, low-level distances and fuel used 
in all mission phases. Such mission planning is very time 
consuming and very susceptible to error since it requires 
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eye-ball interpolation between curves which are as wide as 
the user's pencil lead, and various information is transferred 
from chart to chart. For example, even with extremely precise 
and careful interpretation of take-off distance, an error of 
^ 200 feet could be made due to pencil width alone, exclusive 
of interpolative errors. The combination of these human and 
mechanical errors could build to as much as 500 feet, which 
would be very critical when taking off at maximum gross weight 
on a hot day. Errors in fuel management would be equally as 
critical, if the mission profile dictated return to base or 
ship with minimal fuel. Also, it is noted that utilization 
of the NATOPS performance charts is generally limited to pre- 
flight planning and not feasible for use in high performance 
aircraft, and in several of the helicopters, for' in-flight 
changes in mission profile. A NATOPS pocket check list is 
used for in-flight planning, which is small and therefore 
subject to even higher human and mechanical interpolation 
errors . 

Therefore, a much improved and more accurate method of 
mission planning, both before flight and during flight, is 
needed. A desk top calculator or hand-held computer with 
pre-programmed chips could eliminate these errors during 
pre-flight planning, while in-flight errors could be reduced 
substantially with performance equations programmed into an 
on-board systems computer of the aircraft, or a hand-held 
computer with a dedicated or read only memory (ROM) , a pro- 
gram chip, capable of storing performance equations for all 
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mission phases. Not only would computerization reduce the 
human errors of flight crews who have the NATOPS charts avail- 
able, it would reduce the accident potential of flight crews 
who are airborne without NATOPS and make performance decisions 
from past experience. Such decisions, without accounting for 
all variables which affect aircraft performance, could lead 
to large errors in predicting the capability of the aircraft 
to safely complete a mission phase; and have been cited as 
the cause of a large percentage of Navy /Marine Corps accidents 
in both fixed and rotary wing aircraft. Computerization would 
also greatly enhance shipboard operations by providing in- 
stance Bingo (minimum fuel) profiles to both flight crews and 
the air base, reducing fuel management errors based on exper- 
ience and interpolation of charts. 

This thesis computerizes A-6E aircraft NATOPS take-off and 
cruise performance charts [Ref. 1]. However, the feasibility, 
procedures and techniques are applicable to any aircraft, 
both fixed and rotary wing. The feasibility aspect of the 
thesis is a continuing application of thesis research con- 
ducted for computerization of A-7E performance [Ref. 2]. The 
programs are written in BASIC language for the Hewlett-Packard 
9830A computer, but knowledge of the language is not required 
to demonstrate the feasibility of accurate performance com- 
puterization. The performance equations within the programs 
can be directly utilized in most scientific computer languages, 
with easy adaptation to any computer such as hand-held, desk 
top, or aircraft systems computers. 
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II. EXPERIMENTAL PROCEDURE 



A. GENERAL ORGANIZATION 

Two separate, but equally important, A-6E aircraft per- 
formance charts were evaluated for computerization feasibility. 
The first chart, "Take-off Performance" [Ref. 1, Fig. 11-11], 
was utilized since its family of relatively simple, or low 
order polynomial equations, provided use of basic curves in 
studying and developing the programming techniques. The 
second chart, "Cruise Performance," [Ref. 1, Fig. 11-110], 
provided a more complex family of curves which required more 
extensive analysis in program development. Empirical data 
were taken from the specific NATOPS chart and the Leas t- Squares 
method of polynomial approximation was used to determine the 
coefficients of each curve in the respective family of curves. 

A polynomial expression was then approximated using the same 
method to determine the family of coefficients. In other 
words, polynomial expressions were determined which approxi- 
mated the coefficients of any curve which would fall within 
the bounds of the respective family of curves. This technique 
precluded the need for interpolation between curves. A 
detailed example of the Least Squares method is presented in 
Appendix A, while Figure 1 presents a simplified example of 
approximating a family of curves. 
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Figure 1. Family of curves. 



Least Squares approximations yield polynomial expressions 
for each curve in the form of 
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where for a... • x., i is the specific curve and i is the 

1 j j 1 j r- 

coefficient subscript. 

A polynomial expression for each coefficient family is 
determined in the form of: 



A(0) = • C + b2 • C' 
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A(l) = • C + b 2 • 

A(2) = b^ + b, • C + b- • 

where "C is the curve variable. 

Then, the expression for the curve passing through point 
A in Figure 1 for = C 2 ^ is: 

Y(2.5) = A(o) + A(l) • X 2 5 + A(2) • X 2 5 ^. 

The Least Squares method was programmed for the Hewlett- 
Packard 9830A computer complemented with the Hewlett-Packard 
9266A plotter which was used to verify that the plot of the 
polynomial expressions did, in fact, duplicate the specific 
NATOPS curve. The accuracy of the polynomial approximations 
was determined by evaluating the percentage of deviation of 
the approximated curve at specific points from the same points 
on the actual NATOPS curve. It was possible to achieve greater 
than ninety-seven percent { 91 %') overall accuracy. 

Accuracy of the polynomial expressions for each coefficient 
family was determined by comparing the computed value of the 
coefficient evaluated at each curve of the family with the 
original coefficient approximated for each curve. A more 
detailed description of the computer techniques is contained 
in Ref. [2] . 

B. TAKE-OFF PERFORMANCE 

The A- 6 E aircraft NATOPS normal take-off performance chart 
(Fig. 3] was photographed and enlarged five times its original 
size to allow more accurate measurement of empirical data. 
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The chart was then separated into sub-charts composed of 
families of curves for individual evaluation of computer- 
ization feasibility. The dependent variable or output of 
the first sub-chart became the independent variable or input 
of the succeeding sub-chart. Table I shows the relationship 
of sub-charts to the independent and dependent variables. 



TABLE I 

TAKE-OFF PERFORMAiNCE SUB- CHARTS 



Sub-Chart 


Family of Curves Parameter 


Input 


Output 


I 


Pressure Altitude 


Temperature 


Reference 

Index 


II 


Gross Weight 


Reference 

Index 


No Wind 

Take-off 

Distance 


III 


Take-off Distance 


Wind 


Take-off 

Distance 



Empirical data were taken from the respective sub-chart, 
the coefficient of each curve was determined, and the family 
of curves was plotted to verify accuracy with the original 
sub-chart. Once the points verified for each curve were 
within one percent of the original empirical data, the coef- 
ficients of the family of curves were determined. A computer 
program was written and the coefficients were verified for 
accuracy within two and one-half (2-1/2) to three (3) percent 
of the coefficient previously obtained for each curve. The 
independent variable for the coefficients of the family of 
curves is the specific "family of curves parameter," such as 
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pressure altitude in Sub-Chart I. The dependent variables 
are the coefficients of the curve passing through the specific 
family of curves parameter. These coefficients are then used 
in the output equation of the specific sub-chart. 

In Sub-Chart I, "Pressure Altitude," it was relatively 
simple to determine the coefficients of the family of curves, 
since there were the same linear relationships between each 
curve. But Sub-Charts II and III, "Gross Weight" and "Take- 
off Distance with Wind," had to be broken down into three and 
two sub-families of curves, respectively. This was required 
since the non-linear relationship between curves would pro- 
duce inaccurate coefficients for the family of curves. By 
overlapping smaller sub-families of curves whose relationships 
were more linear, the desired accuracy could be achieved. The 
range covered by each sub-family of curves was determined by 
trial and error iteration until the desired accuracy was 
achieved. Coefficients for the sub-family of curves would 
then be used in the output equation for the specific sub-chart. 
The sub-family coefficients used depend upon the range in 
which the desired family of curve parameter would lie. The 
computer program was written to check for the limits of each 
sub -family range. 

Once satisfactory programs were written for each sub-chart, 
they were integrated into one computer program which would 
yield a final output of take-off distance with wind. The 
final program also includes equations which calculate take-off 
velocity, and unsafe and not-recommended take-off distances. 
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The coefficients for these equations were determined in the 
same manner as individual curves for the sub-charts. A wind 
component computation using basic trigonometry identies was 
used to input winds into Sub-Chart III, and a round-off 
routine was used prior to the final output printout. The 
round-off criterion established for take-off velocity was 
one-half (0.5) knots, while that for take-off distance was 
twenty-five (25) feet. Twenty-five feet was used for take- 
off distance round-off, since at high take-off velocities 
this value would be negligible in rounding down, but would 
still provide a factor of safety for round-up to the nearest 
one hundred (100) feet. 

C. CRUISE PERFORMANCE 

The A-6E NATOPS Cruise Performance Charts are a series of 
separate charts or four phases which require that the output 
of each chart be transferred as the input to the succeeding 
chart until arriving at the final output of "fuel flow." 

These four phases are: 

Phase I Clean Aircraft Transfer Scale 
Phase II Aircraft Reference Number 
Phase III Pounds of Fuel per Nautical Mile 
Phase IV Fuel Flow 

Since Phase I chart consisted of the most non-linear 
curves, in contrast to the linear curves of Phases II, III 
and IV, and in view of project time constraints, only Phase I 
was evaluated for computerization feasibility. This chart 
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was also photographed and enlarged to facilitate extraction 
of accurate empirical data. As in take-off performance, the 
chart was separated into sub-charts for individual evaluation 
of families of curves. This technique varied from that of 
take-off performance in that Sub-Chart II had an overlying 
family of curves Vhich had to be treated as separate entities. 
The dependent variable or output of one curve family became 
an independent variable or input to the second curve family. 
Table II presents the sub-charts with respective curve 
families and variables. 

TABLE II 

PHASE I CRUISE PERFORMANCE SUB -CHARTS 



Sub- Chart 


Family of Curves Parameter 


Input 


Output 


I 


Pressure Altitude 


Gross Weight 


Mach 


II 


Guideline 


Mach 


Transfer 

Scale 


III 


Drag Count 


Transfer 

Scale 


Mach 



Determination of the coefficients of each curve, accuracy 
verification, and determination of curve family coefficients 
were accomplished using the same method described in Take-off 
Performance. Sub-Chart I, "Pressure Altitude," and Sub-Chart 
II, "Guideline," had to be broken down into two overlapping 
sub-families of curves to achieve the desired accuracy. The 
range of each sub-family was determined by trial and error 
until the desired accuracy was achieved. Coefficients for 
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each family of curves were used in the output equation for 
the computer program written for each family of curves. Each 
program with sub-families was checked for the range in which 
the family of curves parameter would lie to ensure that the 
proper coefficients were utilized. 

Since Sub-Chart II had overlying curve families, which 
required that the input variable, Mach number, would intercept 
a guideline at the zero drag count, and proceed down that line 
until intersecting the desired drag count line, a comparing 
and decrementing routine had to be developed. A graphical 
example of this routine is shown in Figure 2. 





Figure 2. Mach comparison and decremation. 
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In this routine, initial Mach (M0) from Sub-Chart I is 
input, as the independent variable, into the zero drag count 
equation which outputs the dependent variable, initial transfer 
scale (T0) . T0 then becomes the independent variable which 
determines the specific guideline to be used by establishing 
its coefficients. The independent variable in the guideline 
equation is Mach, which calculates the dependent variable. 
Transfer Scale (Tl) . T1 is then input as the independent var- 
iable into the selected drag count equation which calculates 
the dependent variable, Mach (Ml). Ml is then compared to 
the previous Mach which was input to the guideline equation. 

For the first iteration. Ml is compared with M0. If the 
absolute value of the difference between the two Machs is not 
J.ess than or equal to 0.01, then M0 is decremented, A = .01, 
and a new transfer scale, Tl ' , is calculated and input back 
into the drag equation to output's new Mach, Ml'. This iter- 
ative and comparison routine is continued until the comparison 
parameter is satisfied. The compared Machs are then averaged 
and output as Best Cruise Mach, M2, while the last computed 
Tl is output as Transfer Scale, T. M2 and T are the final ■ 
outputs of Sub- Chart II. 

Once satisfactory programs were written for each sub-chart 
and verified for accuracy, they were integrated into one 
computer program to yield the final output of Best Cruise Mach 
and Transfer Scale. The numerical values of these outputs 
were not rounded off for aircrew use, since they would be 



19 



transferred as inputs to subsequent phase charts of the cruise 
performance computer program in follow-on thesis work. 
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III. RESULTS AND DISCUSSION 



A. TAKE-OFF PERFORMANCE 

The A-6E NATOPS normal take-off performance chart, Figure 
3, was used to study the feasibility and to expand the basic 
techniques of computerizing NATOPS charts. As previously 
discussed, each sub-chart was evaluated separately and then 
integrated into one computer program. However, Sub-Chart IV, 
Obstacle Clearance, was not utilized in the program since 
this aspect of take-off is seldom critical for the A-6E air- 
craft. Each sub-chart family of curves was relatively simple 
to reproduce using the Least Squares method since the highest 
order of polynomial equation approximated was third order. 

The accuracy of each reproduced curve was verified to be 
greater than ninety-nine percent (99%) . A comparison of Sub- 
Chart II, Gross Weight, of Figure 3 with that plotted from 
the approximated equations. Figure 4, will show that an almost 
exact reproduction of the NATOPS charts can be made by Least 
Squares approximations. The ordinate system of Figure 4 is 
transposed relative to that in Figure 3, since the output 
variable of Sub-Chart I is the input variable to Sub-Chart II. 

While the individual curves of each respective sub-chart 
could be reproduced with a high degree of accuracy, high 
errors were generated when approximating coefficients for the 
family of curves in Sub-Charts II and III. These errors were 
a direct result of the curves in Sub-Chart II becoming more 
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Figure 3. A- 6E Normal Take-off Distance. 
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Figure 4. Take-off Sub-Chart II Approximated Curves. 



non-linear as the gross weight increased, and of the large 
slope changes in Sub-Chart III curves as the take-off dis- 
tances increased. The errors in these sub-charts were 
reduced to less than two percent (2%) by overlapping smaller 
sub-families of curves of more nearly the same linearity and 
slope. The gross weights were reduced to an overlapping 
family of curves reflecting gross weight changes of ten 
thousand (10,000) pounds, while take-off distance was over- 
lapped at the mid-curve of the family. The specific overlap 
points were determined by trial and error until an acceptable 
error tolerance was achieved. 

The final computer program which culminated the respective 
sub-chart programs into one A-6E NATOPS take-off performance 
^^computer program is included in Appendix B. Sample program 
results are in Appendix E, while Table III shows the relative 
accuracy between computer and manual chart determinations of 
take-off distances for representative input parameters. 
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TABLE III 



TAKE-OFF DISTANCE ACCURACY 





1 1 


2 


3 


4 


Gross Weight (lbs.) 


♦ 

34,000 • 


40,000 - 


48,000 - 


53,000 


Pressure Altitude (ft.) 


0 


1,000 


5,000 


0 


Temperature, Deg. (F) 


60 


100 


60 


80 


Runway Heading, Deg. 


360 


360 


360 


360 


Wind Direction, Deg. 


360 


318 


042 


060 


Wind Velocity 


10 


20 


20 


20 


Headwind Velocity 


10 


15 


15 


10 


Take-off Distance: 










Manual Chart 


1,200 


2,550 


4,400 


4,400 


Computer 


1,200 


2,500 


4,300 


4,400 


Accuracy (%) 


100 


98 


97 


100 



B. CRUISE PERFORMANCE 

The A-6E NATOPS Cruise Performance Chart - Phase I is 
shown in Figure 5, while all four phase charts are included 
in Appendix C. The same basic tachniques used in take-off 
performance were applied to cruise performance with exception 
of the iterative routine required for the overlying drag count 
and guideline curves of Sub-Chart II. As in take-off perfor- 
mance, Sub-Chart I, pressure altitude curve family, was rela- 
tively easy to reproduce with individual curve accuracies 
verified as greater than ninety-nine percent (99%) . 
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Figure 5. Cruise Performance, Phase I. 
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However, high errors were generated when approximating coeffi- 
cients for the curve family due to the increased non-linearity 
as pressure altitude increased. These high errors were 
reduced to less than two percent (21) by overlapping sub- 
families of curves from zero (0) to thirty thousand (30,000) 
feet and thirty (30) to forty thousand (40,000) feet. 

In order to achieve accuracies of greater than ninety-nine 
percent (99%) for the individual curves of the drag count 
family in Sub-Chart II, a high number of empirical data points 
were used and the approximated coefficients were limited to 
fourth order. This was required to smooth the curves out and 
prevent oscillations between points. In addition, some points 
in the zero (0) drag count curve were adjusted slightly to 
add additional curve smoothing and better curve reproduction. 
While adjustments were necessary to achieve individual curve 
accuracy, none was required to approximate coefficients for 
the entire family of drag curves, which resulted in greater 
than ninety-eight percent (98%) accuracy. 

Due to the extreme non-linearities of the guideline family 
of Sub-Chart II, the individual curve approximations were 
limited to the range that overlaid only the drag count curves . 
This restriction was required to achieve greater than ninety- 
nine percent (99%) individual curve accuracy, and to limit 
the overlapping sub-family of curves to two. The sub-families 
overlapped at the fourth guideline, as counted from bottom to 
top in Figure 4. 
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Since the Hevvlett- Packard 9830 Computer Plotting Routine 
requires that the X and Y axis increase from the ordinate, 
the transfer scale of Sub-Chart II, Figure 4, had to be in- 
verted to reflect values from zero to thirty, vice thirty to 
zero. This inversion of the Y-axis caused the approximated 
guideline curves to be plotted inverted. Figure 2. However, 
this did not affect the output, transfer scale, since the 
initial point of entry into the guideline curve family is 
only dependent upon the initial transfer scale, which is de- 
rived from the zero drag count equation regardless of curve 
orientation. Also, as discussed previously, the iterative 
routine required that the approximated drag count family of 
curves be transposed. Figure 2, so that Best Cruise Mach could 
be determined by comparing input Mach to the guidelines and 
output Mach from the drag counts. An attempt was made to 
keep the drag count curves oriented the same as the guideline 
curves, but the program would not step out of its iterative 
routine when the input Mach would not pass through the desired 
drag count curve when intercepting the zero drag count curve. 
For example, as depicted in Figure 4, an input Mach of 0.73 
does not pass through a desired drag count of fifty (50). 

The final computer program which culminated the respective 
sub-chart programs into oneA-6E NATOPS Cruise Performance Phase 
I computer program is included in Appendix D. Sample Program 
results are in Appendix E, while Table IV shows the relative 
accuracy between computer and manual chart determinations of 
transfer scale and best cruise Mach for representative input 
parameters . 
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TABLE IV 



PHASE I CRUISE PERFORMANCE ACCURACY 



1 1 


2 1 


3 1 


4 


Pressure Altitude (ft.) 


33,000 


i 

25,000 1 


. 15,000 


37,000 


Gross Weight (lbs.) 


48,000 


50,000 


58,000 


42,000 


Drag Count 


20 


100 


150 


20 


Best Cruise Mach 


i 








Manual 


. 710 


.620 


. 559 


. 718 


Computer 


. 714 


.623 


. 557 


.722 


Accuracy (%) 


99.4 


99.5 


99.6 


99.4 


Transfer Scale 










Manual 


13.75 


9.8 


7.8 


14.4 


i 

j Computer 


13.65 

1 


9.9 


7.9 


14.5 ^ 

1 


■ Accuracy (1) 


99.2 


99. 


98.7 


99.3 i 



IV. CONCLUSIONS AND RECOMMENDATIONS 



Within the scope o£ this thesis project, computerization 
o£ aircra£t Naval Air Training and Operating Procedures 
Standardization Program (NATOPS) £light per£ormance charts 
has been shown to be £easible with a high degree o£ accuracy. 
The computer programs developed £or the A-6E aircra£t are 
adaptable to hand-held calculators, desk calculators, or the 
existing aircra£t systems computer. The £easibility, pro- 
cedures, and techniques discussed are applicable to any air- 
cra£t, £ixed or rotary wing. Utilization o£ computerized 
NATOPS per£ormance charts £or be£ore-£light and during- £light 
mission planning will signi£icantly reduce the human and 
mechanical errors inherent to visual interpretation and inter- 
polation o£ the charts. 

This reduction in chart errors, when coupled with aircrews 
having computer access to make airborne aircra£t per£ormance 
decisions based on all per£ormance variables, rather than on 
past experience only, will greatly enhance sa£ety o£ £light 
in completing any or all mission phases. By enhancing sa£ety 
o£ £light through NATOPS per£ormance chart computerization 
which is accessible to all £light crews, the operational readi- 
ness o£ all Navy/Marine Corps £leet units will be increased. 

Development o£ NATOPS per£ormance charts computer programs 
should be expanded to encompass all phases o£ mission planning 
to include shipboard Bingo (minimum £uel) pro£iles and tactical 
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weapons delivery. These computer programs could be incor- 
porated into hand-held calculators with dedicated read only 
memory (ROM) for use in ready rooms and in aircraft without 
on-board systems computers, such as most helicopters and light 
attack fixed wing aircraft; into desk top computers, capable 
of utilizing cassette tapes to program performance equations 
for all operational Navy/Marine Corps aircraft for use in land 
and ship base operations; and into existing on-board system 
computers, for such aircraft which have the capacity for 
computer software expansion. 

Manual inputs to performance equations incorporated within 
on-board systems computers could be significantly reduced by 
utilizing signals from air data computers for temperature, 
pressure altitude, and velocity inputs; fuel quantity, sensors 
for gross weight inputs; inertial platforms for angle inputs 
such as angle of attack, climb angles and glide paths; and 
weapons selection tapes for drag count inputs. 

Future thesis research should address the interface of 
existing aircraft sub-system output signals as input signals 
to NATOPS performance chart equations incorporated into on- 
board aircraft systems computers. 
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APPENDIX A 



LEAST SQUARES METHOD OF POLYNOMIAL APPROXIMATION 

The Least Squares method of polynomial approximation is 
described in detail in reference [3] . Since data points are 
subject to experimental errors ivhich may be relatively large 
in magnitude, the best curve fit through these points is one 
which minimizes the errors. The Least Squares method of 
approximation minimizes the maximum errors between the data 
point? and the curve by making the sum of the squares of the 
errors a minimum. The errors are measured by the vertical 
and horizontal distances from the points to the line as 
depicted in Figure 6. 




Figure 6. Error Measurement. 

The sum of the squares of the errors is used since the 
function has no derivative at the origin and this could be a 
single large error. The Ler.st Squares method is also in accord 
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with the maximum likelihood principle of statistics. If the 

measurement errors have a normal distribution and if the 

standard deviation is constant for all the data, the sum of 

the squares can be shown to have values of slope and intercept 

which have a maximum likelihood of occurrence. 

t h 

In the development of an n^ degree polynomial which will 
represent the best curve fit through empirical data, n+1 or 
greater pairs of data points are required. Let N represent 
the number of (x,y) data pairs. The functional relationship 

y=a + aTX + a-,x^+....+a x^ (1) 

o 1 Z n 



is assumed with errors defined by: 



e . 

1 



= Y. 
1 



- Y: 



= Y. 
1 



- a 



a, X . 
1 1 



a,,x . 
2 1 



a X . 
n 1 



n 



( 2 ) 



The is used to represent the empirical value correspond- 
ing to x^ of the data pair. The Least Squares criterion 
requires that 



2 ^ 2 ^ 
®1 



'N 



T 2 
Z e . 

i=i ^ 



N 

■E 

i = l 



(Y- - 

o 



a,x. - a-x. 
1 1 2 1 



ax. ) 
n 1 



(3) 



be a minimum. 

The minimum is reached by proper choice of the coefficients, 

a through a . Therefore, the coefficients are the variables 
o ® n 

of the problem. At the minimum for S, the partial derivatives 
of S with respect to the coefficients are zero. Writing the 
equations for the n+1 derivative equation give: 
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0 



N 

3S/8a^ = 0 = Z 2(Y, - - a,x. 

0 0 1 1 



n 



•-a„x.“)(- 1), 



N 



3S/9a, = 0 = Z 2(Y. - a - a.x.- 
1 i-i 1 ° 1 1 



n 



- a^x. )(- xp, 
(4) 



N 

3S/3a„ = 0 = Z 2(Y. - a - a^x. 
n o 1 1 



n 



. . - a X .“ ) ( - X .^ ) 
n 1 ^ ^ 1 ^ 



Dividing each of these equations by -2 and rearranging gives 
the n+1 normal equations to be solved simultaneously: 



'j 

a^N + a, Zx. + a^Zx. +....+ a Zx.^ = ZY., 
oli2i ni i’ 

a zx. + a, Zx. + a^Zx. +....+ a Zx.^'*''^ = Zx.Y., 
oiii Zi ni 11 



2 3 4 

a Zx. +aiZx. +a^Zx. a Zx. 

01 li 2i ni 



n+2 



Z X Y . , 
1 1 ’ 



a Z X + a, Zx.^^^ + a., Z X . , 

0 1 1 1 2 1 



; (5) 

. + a Z X . = Z X Y . , 

n 1 11 



All the summations run from 1 to N. 

Numerical methods are used to solve these equations simul- 
taneously and are discussed in detail in reference [3]. These 
particular equations have an added difficulty in that they have 
the undesirable property known as ''ill-conditioning.” The 
result is that round-off errors in solving them cause unusually 
large errors in the solutions, which are the desired values of 
the coefficients, a^'s. Up to n=4 or 5 the problem is not too 
great, but beyond this point special orthogonal polynomial 
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methods are needed, which is a separate topic in reference [3]. 
However, from the experimental point of view, functions more 
complex than fourth- degree are rarely needed, and when they 
are, the problem can often be handled by fitting a series of 
polynomials to subsets of data. 

To illustrate the use of equation (^5), a quadratic is fit 
to the data of Table V. 



TABLE V 
CURVE FIT DATA 

Data Pairs 



X . 
1 


0.0 


0.2 


0.4 


0.7 


0.9 


1.0 


Y. 

1 


1.016 


0.768 


0.648 


0.401 


0.272 


0.193 



Summations for N=6 



X . 
1 


2 

X . 

1 


X . 

1 


X.4 

1 


Y. 

1 


x.Y. 
1 1 


2 V 

X . Y . 
1 1 


3.2 


2.5 


2.144 


1.9234 


3.298 


1.1313 


0.74421 



The data are actually a perturbation of the relationship 
y=l-x+0.2x^ (6) 

The curve of equation (6) and the perturbated data are 
shown in Figure 7. 
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Figure 7. Curve fit data plot. 

The summations tabulated in Table V are substituted into 
equations (5) to give the following set of equations: 

6 a^ + 3.2 a^^ + 2,5 = 3.298, 

3.2 a^ + 2.5 + 2,144 a^ = 1.1313, (7) 

2.5 a^ + 2.144 a^ + 1.9234 a2 = 0.74421. 

Using numerical methods of solution to simultaneous 
equations yield: 

a = 0.9986, a. = -1.0060, a. = 0.2103. 
o ’1 ’2 

Therefore, the equation of the Least Squares approximated 
polynomial through the data of Table V is: 

y = 0.999 - 1.006 x + 0.210 (8) 
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The coefficient of equation (8) is slightly different from 
that of equation (6) because of the errors in the empirical 
data . 
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APPENDIX B 



TAKE-OFF PERFORMANCE COMPUTER PROGRAM 
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APPENDIX C 



CRUISE PERFORMANCE PHASE CHARTS 

NAVAIR 01-85ADF-1 



CRUISE PERFORMHCt 



PERFORMANCE 
Mission Planning 




PHASE ! - CLEAN AIRCRAFT TRANSFER SCALE 



DATE 15 MARCH 1971 
DATA BASIS; ESTIMATED 



FUEL GRADE; JP-5 

FUEL DENSITY; 6.8 LB./GAL. 




AOf]-311 



Figure 11-11 

\l 



PERFORMANCE 
Mission Planning 



NAVAIR 01-85ADF-1 



CIWISC PERFORmm 





PHASE II ^ AIRCRAFT REFERENCE NUMBER 



DATE 15 MARCH 1971 
DATA BASIS: ESTIMATED 



FUEL GRADE: JP.5 

FUEL DENSITY; 6.8 LB./GAL 




ADfl-012 



Figure 11-12 



11-18 
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NAVAIR 01-85ADF-1 




CRUISC PERfORfmCC 




PERFORMANCE 
Mission Planning 




PHASE III - POUNDS OF FUEL PER NAUTICAL MILE 



DATE 15 MARCH 1971 
DATA BASIS: ESTIMATED 



FUEL GRADE: JP-5 

FUEL DENSITY: 6.8 LB./GAL. 












AD#1-313 



Figure 11-13 



11-19 
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PERFORMANCE 
Mission Planning 



NAVAIR 01-85ADF-1 



1 



PHASE IV - FUEL FLOW 



CRWSf PElifORMAHa 





DATE 15 MARCH 1971 
DATA BASIS: ESTIMATED 



FUEL GRADE: JP*5 

FUEL DENSITY: 6.8 L8./GAL. 




A0ft4U 



Figure 11-14 



11-20 
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APPENDIX D 

CRUISE PERFORMANCE COMPUTER PROGRAM 
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APPENDIX E 



SAMPLE PROGRAM RESULTS 
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PRESSURE RLTITUDE = 33000 

GROSS NEIGHT -- 48000 

DRRG COUNT = 20 

C R U I S E f’1 f=l C H = 0 . 7 1 3 8 5 9 6 5 0 

T R R N S F E R S C R L E = 13, 3 4 g 0 2 0 0 5 
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